Competition between spin and charge polarized states 
in nanographene ribbons with zigzag edges 
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Effects of the nearest neighbor Coulomb interaction on nanographene ribbons with zigzag edges are 
investigated using the extended Hubbard model within the unrestricted Hartree-Fock approxima- 
tion. The nearest Coulomb interaction stabilizes a novel electronic state with the opposite electric 
charges separated and localized along both edges, resulting in a finite electric dipole moment pointing 
from one edge to the other. This charge-polarized state competes with the peculiar spin-polarized 
state caused by the on-site Coulomb interaction and is stabilized by an external electric field. 

DOI:10.1103/PhysRevB PACS number(s):71.10.Hf, 73.22.-f,73.20.At,77.22.Ej 



New materials composed of carbon atoms have been 
attracting much attention both in the fundamental sci- 
ence and in the potential application to nanotechnology 
devices [1-3]. One of the main interests is the variety 
of their physical and chemical natures depending on ge- 
ometries [1-9]. In carbon nanotubes, diameters and chi- 
ral arrangements of hexagonal pattern on tubules decide 
whether they are metallic or not [1-3] . Besides the closed 
carbon molecules like carbon nanotubes, carbon systems 
with open boundaries also show unusual features. 

In nanographencs, i.e., nanoscale graphites, their pe- 
ripheral edges strongly affect the electronic states [4]: A 
nanographene ribbon with zigzag edges (zigzag ribbon) 
has the nonbonding molecular orbitals localized mainly 
along the zigzag edges (edge states). They constitute 
partly flat bands and cause a sharp peak in the density of 
states at the Fermi energy. On the other hand, one with 
armchair edges does not have such edge states. When the 
on-site Coulomb interaction is introduced within the un- 
restricted Hartree-Fock (UHF) approximation, it causes 
the Fermi-instability of the flat bands and forms the spin- 
polarized (SP) states: most spin density concentrates 
along the zigzag edges forming the ferrimagnetic orders 
due to the edge states. As to the effects of the Coulomb 
interactions, there still remains the qualitatively impor- 
tant problem of the competition between the spin and 
charge orderings in this system due to the on-site and 
nearest neighbor Coulomb interactions, respectively. In 
fact, for two dimensional (2D) boundless graphite, the 
possibility of the stable charge density wave (CDW) state 
compared with the spin density wave (SDW) state was 
reported [10,11]. 

In this paper, we demonstrate that the nearest 
neighbor Coulomb interaction stabilizes a novel charge- 
polarized (CP) state with a finite electric dipole moment 
in zigzag ribbons and it competes with the SP state. 

Figure 1 illustrates the geometry of a zigzag ribbon 



with an inversion symmetry. N and L are the width and 
length of the ribbon, respectively. The dashed rectangle 
denotes a unit cell. The periodic boundary condition is 
set along the y axis. Since the zigzag ribbon is bipartite, 
A and B sites are assigned by black and white circles, 
respectively. Note that all the twofold coordinated sites 
in the lower and upper zigzag edges belong to the A and 
B sublattices regardless of N, respectively, at which edge 
states are mainly localized. 

We treat a half-filled it electron system on the zigzag 
ribbon with the extended Hubbard Hamiltonian, 

Heh = -t (4,s^s + H.c.) 

1 (id) 

Here denotes a pair of neighboring carbon sites and 
—t is the transfer integral between them. The operator 
cj s creates a tt electron with spin s on i-th site and Ci yS 
annihilates the one. n^ s is the corresponding number 
operator and ra = X) s n M- The parameters U and V 
denote the on-site and nearest neighbor Coulomb inter- 
actions, respectively. 

We adopt the standard UHF approximation [4] and the 
mean field Hamiltonian is 

h = - ]T [(< + y(4 sCl , s ))4 sCj , s + B.c.-y|(4 sCj , s )| 2 ] 
+ u E<E ( Ui '- s ~ 2> ni ' s ~ ( n *d)( n a) + j ) 

is ~ 

+ Vj2(( n i - l ) n o + ( n J ~ l ) n i - ( n i)( n o) + !)• ( 2 ) 

(id) 

The charge density at i-th site is given by d{ = 1 — (ri.;) 
while the z component of spin density at i-th site is given 
by s, = ((rii >T ) - (n a ))/2. 
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Figure 2 (a) shows the charge density distribution in 
the CP state with NxL = 4x20, where U = OM and 

V = OAt. For comparison, Fig. 2 (b) shows the spin 
density profile of the SP state for U — t and V = 0. 
The CP state has no spin density at every site, while the 
SP state has no charge density at every site. In the CP 
state, the lower zigzag edge are charged positively, while 
the upper zigzag edge are charged negatively. The dis- 
tribution pattern of the charge density in the CP state 
is just like that of the spin density in the SP state. The 
charge polarization in the CP state is explained by the in- 
terplay between the edge states and the Fermi-instability 
of the flat bands due to V as follows, just as the spin po- 
larization explained by the edge states and U [4]. First, 

V stabilize the charge ordered state similar to the CDW 
state with A and B sites charged positively and nega- 
tively, respectively, that is the Fcrmi-instability of the 
flat bands due to V. Therefore, all the twofold sites in 
the lower zigzag edge, A sites, are charged positively, 
while those in the upper zigzag edge, B sites, are charged 
negatively. Moreover, magnitudes of the charge densities 
at twofold sites are much larger than those at adjacent 
threefold sites with unlike charges due to the edge states. 
Therefore, both zigzag edges are charged oppositely and 
the CP state is regarded as a peculiar CDW state caused 
by the edge states. 

The CP state has a finite electric dipole moment point- 
ing from one edge to the other along the x-axis. It is 
formed in the homogeneous system and breaks the in- 
version symmetry. The x component of polarization op- 
erator is P x = ea^ ji Xi(l — n*), where e and a are the 
elementary electric charge and the C-C bond length in 
a zigzag ribbon, respectively. Xi is the x coordinate of 
the i-th site in unit of a. The electric dipole moment per 
unit cell \i = (P x )/(L/2) in Fig. 2 (a) is 0.88 [ea] (= 5.9 
Debye for a = 1.4 A). 

Figure 3 shows the phase diagram in the parameter 
space of U and V for the competition between the CP 
and SP states, where NxL = 4x40. Above the phase 
boundary the CP state is lower in energy than the SP 
state, otherwise the SP state is lower. The boundary is 
determined as the crossing points of their energy curves 
by the numerical interpolation. Near the boundary, the 
CP and SP states coexist for given parameter sets, show- 
ing the first order phase transition. 

At U ~0, the boundary rises up with the infinite gradi- 
ent, indicating that the CP state is higher in energy than 
the SP state in the weak coupling limit. Since Fig. 2 
suggests that the magnitudes of charge densities become 
vanishingly small at threefold sites while they remain fi- 
nite at twofold sites, the CP state cannot be stabilized by 

V and is simply destabilized by U. This comes from the 
edge states absent in nanographene ribbons with arm- 
chair edges and in 2D boundless graphites. 

When V exceeds V c ~t/2 below the phase boundary, 
the CP state gradually loses their localized nature and 



the magnitudes of charge densities grow continuously at 
threefold sites. In the strong coupling limit, the mag- 
nitude of the charge density becomes almost the same 
among all sites and the CP state turns into the CDW 
state, showing the crossover from the CP state to the 
CDW state. This is just like the crossover from the SP 
state to the SDW state at U c ~t above the phase bound- 
ary [4]. In the CDW state, the magnitudes of unlike 
charge densities are the same so that the total charge 
along each edge cancels out, respectively, while (P x ) re- 
mains finite due to the non-uniform width of the zigzag 
ribbons. In the continuous change to the CDW state, 
(P x ) increases to approach LN/A. 

The inset of Fig. 3 shows the phase diagram in a 
reduced scale. In the strong coupling limit, the phase 
boundary between the CDW and SDW states approaches 
the line U — ijV derived from equalizing their Coulomb 
energies per site Eqdw = U/4— Vrj/2 and Esdw = 
—U/A, where r] = 3 — 1/N is the average number of adja- 
cent sites which originates from threefold sites in the unit 
cell with 27V sites except for two twofold sites at edges. 
Note that the CDW-SDW phase boundary in ID or 2D 
system is generally close to the line with r\ being the 
number of adjacent sites for each system [11-17], while 
the edge states in a zigzag ribbon cause a new feature, 
i.e., the deviation of the phase boundary from it in the 
weak coupling limit. In the wide ribbon limit, the phase 
boundary approaches U = 3V for 2D graphite since the 
weights of the edges states diminish and 77 approaches 3. 

Finally, we demonstrate how an external electric field 
stabilizes the CP state. Since the external field couples 
with /z, the CP state with finite electric dipole moment 
/icp is effectively stabilized while the SP state with van- 
ishing electric dipole moment /isp not. Under the static 
electric field with the amplitude E ext polarized along the 
x axis, Hamiltonian is given by H = Heh — E ext P x . 
The unit of E ext is t/(ea) which is 7.1xl0 7 V/cm for 
t = 1 eV. Figure 4 (a) shows the total energies of the 
CP and SP states as functions of E ext for U = 0.9t and 
V = 0.6t. At E ext — 0, the SP state is lower in energy 
than the CP state. Applying the electric field, the SP's 
energy depends parabolically on E ext with the peak at 
the origin and scarcely decreases. This is because the ini- 
tial vanishing ^sp eliminates the first order perturbation 
energy. Concomitantly, the induced fxsp is proportional 
to E ext as shown in Fig. 4 (b). On the other hand, the 
CP's energy depends linearly on E ext and decreases ap- 
preciably, because the initial fic p is much larger than the 
induced one and the first order energy correction is dom- 
inant. Note that the CP state becomes lower in energy 
than the SP state when E ext >A.5xlO~ 3 [t/(ea)}. Figure 
4 (b) shows discontinuous change of \x at the transition, 
due to the first order one. 

The followings are devoted to discussions. First, the 
nanographene ribbon is regarded as a quasi ID system 
where quantum fluctuations are important and an UHF 
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solution is not a reliable ground state. Such a system gen- 
erally has several UHF solutions with almost degenerate 
energies and the natures of their long range orders appear 
as short range orders in the correlated ground state [18]. 
On the other hand, the parameter set U = 2.5t and V = t 
for 2D graphite borrowed from the coronen molecule [11] 
falls into the SP side in the vicinity of the phase bound- 
ary in Fig. 3. This suggests that the CP and SP states in 
zigzag ribbons are close in energy and there will be the 
charge and spin fluctuations along zigzag edges reflecting 
the CP and SP states. In fact, Hikihara et al [5] showed 
that the correlated ground state in the Hubbard model 
for zigzag ribbons is the spin-singlet state without sym- 
metry breaking while the ferrimagnetic correlation at the 
zigzag edges remains robustly. We expect that the CP 
state is a stable UHF solution as the SP state is proved to 
be in the random phase approximation [6] . The symme- 
try breaking /icp is due to the UHF approximation and 
the inversion symmetry will be restored in the correlated 
ground state with fi = 0. It will include the degener- 
ate CP states with oppositely signed fj,cp as important 
configurations and exhibit large quantum fluctuation of 
fi similar to the quantum paraelectricity [19,20]. 

Secondly, the electron correlation will change the na- 
ture of the phase transition obtained in the UHF approx- 
imation. In a 2D extended Hubbard model with square 
lattice, the quantum Monte Carlo study showed that the 
transition is similar to that in ID system [12], where ex- 
tensive studies clarified that in the weak and strong cou- 
pling limits there are continuous (second order) and dis- 
continuous (first order) phase transitions at U ~ 2V, re- 
spectively [13-17]. This may hold for the quasi ID zigzag 
ribbon, while the divergent gradient of the boundary at 
U~0 will survive, i.e., the charge fluctuations reflecting 
the CP states will be suppressed by the edge states in 
the weak coupling limit under the correlation since the 
edge states are robust against the correlation [6,5]. 

Thirdly, this model is based on the screened Coulomb 
interactions that effectively take into account the corre- 
lation beyond the UHF approximation as seen in studies 
of polyacetylene [21]. As to the effect of the long range 
Coulomb interactions, the ah initio study in 2D graphite 
suggests that it prefers the CDW state [10]. However, 
it may destabilize the CP state since the long range 
Coulomb repulsions between like charges along the same 
edges prevail those attractions between unlike charges 
along the opposite edges. If the excess repulsive energy 
surpasses the energy gain due to the charge ordering in 
this model, the CP state will be destabilized. Even if it is 
the case, the CP state is stabilized under the electric field 
as shown in Fig. 4 and will affect the dielectric response 
of the system. The effects of the long range Coulomb 
interactions and the correlation will be studied in future. 

Fourthly, the effect of the geometrical topology on the 
electronic state is another interesting problem. Recently, 
Tanda et al. synthesized Mobius strips formed by the 



single crystals of NbSe3 [22]. This work urged theo- 
retical studies on the effects of the Mobius boundary 
condition [7,23,24]. In the zigzag ribbon, this topology 
causes a ferrimagnetic domain wall between the up and 
down spin domains localized along zigzag edges within 
the Hubbard-UHF model [7]. V induces another novel 
state in a Mobius strip: the ferroelectric domain wall 
connecting the domains with the opposite electric dipole 
moments [25]. The competition between these states in 
the exotic topology will be studied elsewhere. 

Finally, the present results are not only applied to the 
nanographene ribbons but also to boron-nitride nano- 
ribbons [26,27]. We conjecture that the novel CP state is 
valid for the general 2D tt electron systems with zigzag 
edges and affects their dielectric response to be different 
from those without zigzag edges. 

In summary, we found the novel charge-polarized state 
with a finite electric dipole moment caused by the in- 
terplay between the nearest neighbor Coulomb interac- 
tion and the edge states in nanographene ribbons with 
zigzag edges using the extendcd-Hubbard-UHF model. 
This state competes with the spin-polarized state caused 
by the on-site Coulomb interaction and the former is ef- 
fectively stabilized by an external electric field. 

This work has been supported partly by Special Co- 
ordination Funds for Promoting Science and Technology 
and by NEDO under the Nanotechnology Materials Pro- 
gram. 

FIG. 1. Schematic structure of a bipartite nanographene 
ribbon with zigzag edges (zigzag ribbon); •, A site; o, B site. 
Dashed rectangle denotes a unit cell. 



FIG. 2. (a) Charge density distribution, di, of the 
charge-polarized (CP) state for U = 0.3t and V = OAt, and 
(b) the z component of spin density distribution, s,, of the 
spin-polarized (SP) state for U = t and V = 0, on a zigzag 
ribbon with 4x20 sites; •, positive density; o, negative den- 
sity. Radii of circles denote the magnitudes of the charge or 
spin densities; the maximal one is 0.18 in (a) and 0.12 in (b). 



FIG. 3. Phase diagram in the parameter space of U and V 
for a zigzag ribbon composed of 4x40 carbon atoms. Solid 
curve is the phase boundary between the CP and SP states 
which fits numerical data (•). The inset shows the phase 
diagram in a reduced scale. The dashed line denotes the 
CDW-SDW phase boundary in the strong coupling limit. 



FIG. 4. External static electric field E ext dependence of (a) 
energies, and (b) electric dipole moment fj, of the CP state (•) 
and SP state (o). 
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